Negative Magnetoresistance of Granular Metals in a Strong Magnetic Field 
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The magnetoresistance of a granular superconductor in a strong magnetic field destroying the 
gap in each grain is considered. It is assumed that the tunneling between grains is sufficiently large 
such that all conventional effects of localization can be neglected. A non-trivial sensitivity to the 
magnetic field comes from superconducting fluctuations leading to the formation of virtual Cooper 
pairs and reducing the density of states. At low temperature, the pairs do not contribute to the 
macroscopic transport but their existence can drastically reduce the conductivity. Growing the 
magnetic field one destroys the fluctuations, which improves the metallic properties and leads to the 
negative magnetoresistance. 

PACS numbers: 73.23.-b, 74.80.Bj, 74.40. +k, 72.I5.Rn 



In a recent experiment jjj, the resistance of a system 
of Al grains as a function of the magnetic field was stud- 
ied. The samples were quite homogeneous with a typical 
diameter of the grains 120 ± 20A and the grains formed 
a 3-dimensional array. Destroying the superconductivity 
by the magnetic field the authors could observe that in 
a sufficiently strong magnetic field the system had a fi- 
nite resistivity. The applied magnetic fields reached 17T, 
which was more than sufficient to destroy also the super- 
conducting gap in each grain. 

The dependence of the resistivity on the magnetic field 
observed in Ref. |Q was not simple. At strong magnetic 
fields, the resistivity increased when decreasing the mag- 
netic field. Only at sufficiently weak magnetic fields the 
resistivity decreased and finally the samples displayed su- 
perconducting properties. A similar behavior had been 
reported in a number of publications |^,^|. A negative 
magnetoresistance due to weak localization effects is not 
unusual in disordered metals 0]. However, the magne- 
toresistance of the granulated materials is quite notice- 
able in magnetic fields exceeding 10T, which is many or- 
ders of magnitude higher than the typical values relevant 
for the weak localization. 

The aim of this Letter is to demonstrate that the mag- 
netoresistance of a granulated metal in a strong mag- 
netic field and at low temperature must be negative due 
to superconducting fluctuations. We consider a system 
of granules coupled to each other. The tunnelling ampli- 
tude between the grains is assumed to be large, such that 
the system without electron-electron interactions would 
be macroscopically a good metal and would not be sen- 
sitive to a magnetic field. The superconducting gap in 
each granule is assumed to be suppressed by the strong 
magnetic field. All the interesting behavior considered 
below originates from the superconducting fluctuations 
that lead to a suppression of the density of states (DOS) 
but do not help to carry an electric current. 

Theory of superconducting fluctuations near the tran- 



sition into the superconducting state has been developed 
long ago [||-f7j (for a review see Ref. ||). Above the 
transition temperature T c , non-equilibrium Cooper pairs 
are formed and a new channel of charge transfer opens 
0. Another fluctuation contribution comes from a co- 
herent scattering of the electrons forming a Cooper pair 
on impurities ||. Both the fluctuation corrections in- 
crease the conductivity and lead to a positive magne- 
toresistance. Formation of the non-equilibrium Cooper 
pairs results also in a fluctuational gap in the one-electron 
spectrum Rj but in conventional superconductors the 
first two mechanisms are more important and the con- 
ductivity increases when approaching the transition. A 
small decrease of the transverse conductivity is possible 
in layered materials |J in a temperature interval not very 
close to the transition. It is relevant to emphasize that 
the study of the fluctuations has been done near the crit- 
ical temperature T c in a zero or a weak magnetic field. 

In granulated materials, the superconducting gap in 
each granule can be destroyed at low temperature by a 
strong magnetic field. At magnetic fields H exceeding the 
critical field H c virtual Cooper pairs can still be formed. 
However, as it will be shown below, the influence of these 
pairs on the macroscopical transport is drastically differ- 
ent from that near T c . The existence of the virtual pairs 
leads to a reduction of the DOS but, in the limit T — > 0, 
these pairs cannot travel from one granule to another. 
As a result, the conductivity a can be at H > H c con- 
siderably lower than conductivity do of the normal metal 
without an electron-electron interaction. It approaches 
the value do only in the limit H 3> H c , when all the su- 
perconducting fluctuations are completely suppressed by 
the magnetic field. 

For explicit calculations we assume that the super- 
conducting pairing inside the grains is mainly destroyed 
by the orbital mechanism. The Zeeman splitting is not 
important provided the radius R of the grains is large 
enough, so that R ^> £ (eqt)^ 1 , where £ is the supercon- 
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ducting coherence length, eg is the Fermi energy and r 
is the elastic mean free time. This condition is well sat- 
isfied in grains with R ~ 100 A studied in Ref. |Q. This 
limit is opposite to the one considered recently in Ref. 

where the Zeeman splitting was assumed to be the 
main mechanism of destruction of the Cooper pairs. 

We assume that electrons can hop from grain to grain 
and the tunneling energy t is in the interval 

S < t < A (1) 

where 6 — (vV) 1 is the mean level spacing in a single 
granule, v = mpo/2ir 2 is the DOS of the metal in the ab- 
sence of interactions, V is the volume of the granule, and 
Ao is the BCS gap at T = in the absence of a magnetic 
field. The left inequality in Eq. (|l|) guarantees that the 
system is macroscopically a good metal jllj and localiza- 
tion effects can be neglected. Moreover, charging effects 
are also not important in this limit. In the limit R <C £ 
considered below the superconducting fluctuations in a 
single grain are effectively zero-dimensional (0D). 
The Hamiltonian H of the system can be written as 

H = H Q +H T , (2) 

where Hq is a conventional Hamiltonian with an electron- 
phonon interaction in the presence of a strong magnetic 
field. The term Ht in Eq.(|^) describes the tunneling 
from grain to grain and has the form (see e.g. Ref. [fl2"||) 

Ht = Ujp g alp a 3Q + h - c - ( 3 ) 

where a\ p {a,i p ) are creation (annihilation) operators for 
an electron the grain i and state p. 

Correspondingly the a.c. current density j (i) is writ- 
ten using linear response formulae as 
t 

j(t)=i J <[j(0,i),j(0,i')])A(O< (4) 

— oo 

where j is the tunneling current operator and the angle 
brackets stand for averaging over both quantum states 
and impurities in the grains. In principle, the grains can 
be clean and the electrons can scatter mainly on the sur- 
face of the grains. However, provided the shape of the 
grains corresponds to a classically chaotic motion of the 
electrons, the clean limit should be described in the 0D 
case by the same formulae. 

We carry out the calculation of the conductivity mak- 
ing expansion both in fluctuation modes and in the tun- 
neling term Ht- The right part of Eq. (Q) can provide a 
small parameter of the expansion in the tunneling. As in 
conventional superconductors, we can single out contri- 
butions due to corrections to the DOS and Aslamazov- 
Larkin (AL) oal and Maki-Thompson (MT) omt correc- 
tions. Diagrams describing these contributions are rep- 
resented in Fig. 1. As a result, the total conductivity a 
can be written as 



o = (TDOS + a AL + <?MT (5) 
where udos is given by equation 

+OC 

<tdos=<Jo(4T)- 1 f [p(e)/u ] 2 co8h- 2 (^=)de, (6) 



(To = 4-/Te 2 <i~ 1 (t/S) is the classical conductivity of the 
granular metal and v (e) is the density of states modified 
by the superconducting fluctuations and averaged over 
the impurities. 

The main correction 8v (e) to the DOS of the non- 
interacting electrons v§ is described by the diagram in 
the Fig. la, while the terms (Jmt and oal are given 
by Figs, lb and lc, respectively. The calculation of the 
diagrams can bee performed for the Matsubara frequen- 
cies e n — 7rT(2n+ 1) using temperature Green functions. 
At the end one should, as usual jl3| , make the analytical 
continuation ie n — » e . The magnetic field can be consid- 
ered in the quasiclassical approximation, which results in 
additional phases in Green functions. 

The diagrams in Fig.l contain essentially the averaged 
one-particle Green functions, the impurity vertices C and 
the propagator of the superconducting fluctuations K. 
The functions C and K depend on the coordinates and 
time slower than the averaged one-particle Green func- 
tions. As a result, the magnetic field affects only the 
vertex C and the propagator K, whereas the phases of 
the Green functions drawn in Fig.l outside these blocks 
cancel. So, reading the diagrams in Fig.l one should re- 
place the solid lines by the functions 

(ie n , p) = (ie n - £ (p) + i (2T)- 1 sgne^j ~* (7) 

The impurity vertex entering these diagrams is equal 
to {2-ki>ot) 1 C (ie n) iflk — i£ n ), where C is the so-called 
Cooperon. It obeys the following equation 

(d HV+ (2e/c) A) 2 + |2e„ -Sl k \)c (r, r') (8) 
= 2tw <5 (r - r') 

where Dq = Uqt/3 is the classical diffusion coefficient. 
The vector-potential A (r) should be chosen in the Lon- 
don gauge. If the shape of the grain is close to spherical, 
the vector-potential is expressed through the magnetic 
field H as A (r) = [H x r]/2. 

All relevant energies in the problem are assumed to 
be much smaller than the energy of the first harmonics 
E c = Dqtt 2 /R 2 playing the role of the Thouless energy 
of a single grain and this allows to keep only the zero 
harmonics in the spectral expansion. One can find the 
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eigenvalue £o (H) of this harmonics using the first order 
of the standard perturbation theory 



So (H) = (2e/cf D < A 2 > 



(9) 



where < ... >o stands for the averaging over the volume 
of the grain. For a grain of a nearly spherical form one 
obtains 



5 \ir<po 



(10) 



where tf>o — 7rc/e is the flux quantum and <f> is the mag- 
netic flux through the granule. 

Within the zero-harmonics approximation, the func- 
tion C does not depend on coordinates and equals 

C{ie ni iVL k -ie n ) =2^u () {\2e n ~m k \+£o{H))- 1 (11) 

The propagator of the superconducting fluctuations is 
calculated summing, as usual, the superconducting lad- 
der containing the products Gi £n Gi^ k _j En . Neglecting 
weak localization corrections one reduces the averaging 
of the propagator K over the impurities to the averag- 
ing of these products that give finally C . The tunneling 
can be neglected here provided the right part of Eq. (|l|) 
is fulfilled. However, the tunneling from grain to grain 
should be taken into account when calculating K. This 
can be done making expansion in H T , Eq. (§). Each now 
vertex arising in this expansion should be also dressed by 
impurity lines. 

Although the final result can be written for arbitrary 
T and H, let us concentrate on the most interesting case 
T <C T c , H > H c . Assuming for simplicity that the 
granules are packed into a cubic lattice and using the 
momentum representation with respect to coordinates of 
the grains we obtain in this limit 



K(iD, k ,q) = -v' Q 



In 



£o(H) 
A 



£o(H) 



■»?(q) 



7?(q) = (4/3tt) J iS/So(H)) (1 - cos qd) ■ 



(12) 



where J = (tt 2 /4) (t/S) i , q is the quasi-momentum, 
d = 2R, and |f2fc| <C £q{H). The pole of the propa- 
gator K, Eq. (p"2|), at q = 0, Q k = determines the field 
H c , at which the BCS gap disappears 



£ {H C ) = A 



(13) 



The result for iJ c , Eqs. ( ]To| , p"3| ), agrees with the one 
obtained long ago by another method The term 

r)(q) in Eq. (|l2|) is very important if H is close to H c . For 
e n > the contribution of the diagram Fig. la, is: 



5v(ie n ) 



2iT 
vo 



J2 / K{in kl ci)C 2 ( 



d 3 q 
(14) 



After calculation of the sum over Q k in Eq. (|l4|), one 
should make the analytical continuation is n — > e. At low 
temperatures it is sufficient to find 5v = 5v (0). 
Remarkably, Eqs. 
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do not contain explicitly 
the mean free time r. This is a consequence of the zero- 
harmonics approximation, which is equivalent to using 
the random matrix theory (RMT) plj . (The parameter 
r enters only Eq. ( |l0| ) giving the standard combination 
£ o (H) describing in RMT the crossover from the orthog- 
onal to the unitary ensemble). This justifies the claim 
that the results can be used also for clean grains with a 
shape providing a chaotic electron motion. 

Using Eqs. (^, |l2}|l4| ) one can easily obtain an ex- 
plicit expression for <jdos for H — H c <C H c . In this 
limit, one expands the logarithm in the denominator of 
Eq.(|l2"|) and neglects the dependence of C on e n and 
fife because the main contribution in the sum over fl k 
comes from £l k ~ £o(H) — £q(H c ) <C Ao. The result for 
8&dos = °o — &DOS can be written as 



^0 



25 f -7T- 1 

a7\£ 



< lnr?(q) > q 



< V \q) > 



T/A « fj 
< T/A < 1 



(15) 



»j(q) = f?(q) + 2/i, 



< ••• > c 



V / (...)dq/(27r) 3 



where h = (H — H c ) / H c . We see that the correction 
to the conductivity is negative and its absolute value de- 
creases when the magnetic field increases. The correction 
reaches its maximum at H — > H c . At zero temperature, 
the maximum value of Sodos/vq is of order S/Aq. As 
temperature grows, the correction can become larger and 
reach for T ~ Ao the order of magnitude of J -1 . Both 
the values are smaller than unity because we work in 
the region specified by Eq. (|l|) and this justifies the dia- 
grammatic expansion we use. The correction can become 
comparable with ao when J ~ 1 , which would mean that 
we are not far from the metal-insulator transition. For 
such values of J one can use Eq. (|iq ) only for rough 
estimates. Apparently, the parameters of the samples of 
Ref. correspond to the region J ~ 1 ; <5/A ~ 1. 

In the opposite limit H H c the correction to 
(To can still be noticeable. In this case we can ne- 
glect the dependence of the superconducting propagator, 
Eq. (|l2[) , on il k and on the tunneling term (because now 
ln(£ Q (H)/£ (H c )) > 1) and finally obtain 

Sctdos/cto - -(1/3) (S/£o(H)) In" 1 (£ (H)/A ) , (16) 

which means that in the region H ^> H c the correction 
to the conductivity decays essentially as Saoos ~ H~ 2 . 
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Let us emphasize that the correction to the conduc- 
tivity coming from the DOS remains finite in the limit 
T — > 0, thus indicating the existence of the virtual 
Cooper pairs even at T — 0. 




FIG. 1. Diagrams describing corrections to conductivity 
due to superconducting fluctuations. The wavy lines denote 
the propagator of the fluctuations, dashed lines describe im- 
purity scattering. 

In order to calculate the entire conductivity, Eq. (||), 
we must investigate the AL and MT contributions (Figs, 
lc and lb). In conventional superconductors near T c , 
these contributions are most important leading to an in- 
crease of the conductivity. 

In the granular materials the situation is much more 
interesting. It turns out that both the AL and MT con- 
tributions vanish in the limit T — > at all H > H c and 
thus, the correction to the conductivity comes from the 
DOS only. Leaving the presentation of the details for a 
future publication [p| we list here only the final results. 

The AL correction to the conductivity gal corre- 
sponds to the tunneling of the virtual Cooper pairs and 
must be proportional to t in contrast to the one-electron 
tunneling determining gdos- All calculations are done 
in the same approximation as when calculating gdos- ln 
the limit TcT c , H — H c <C H c , the result for gal can 
be written as 



gal 47T 2 i 2 T 



<A(q)[7j(q)] WfttZ >, (17) 



00 



9Ag 



where A (q) = 4ttT (3A ^(q)) for T < A ?) and 
A (q) = 1 for T > A ?7 . 

As usual, the MT diagrams have both regular and 
anomalous part. For the problem considered, they are 
of the same order of magnitude but have opposite signs. 
Moreover, at T = they cancel each other. The final 
result for the MT contribution takes the form 



gmt 8tt 2 T 2 6 



00 



9Ag 
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E 

i=l 



< B (q) cos qid > q 



(18) 



where B (q) = — ir 1 ln?)(q) for T <C Aqtj and B (q) = 
2T(A 7y(q))- 1 forT> A 7y. 

The temperature and magnetic field dependence of 
gal and gmt is rather complicated but they are defi- 



nitely positive. The competition between these correc- 
tions and gdos determines the sign of the magnetoresis- 
tance. We see from Eqs. ( fL7[ |l8| ) that both the AL and 
MT contributions are proportional at low temperatures 
to T 2 . Therefore the gdos in this limit is larger and the 
magnetoresistance is negative for all H c . In contrast, at 
T ~ T c the AL and MT corrections can become close to 
H c larger than gdos resulting in a positive magnetoresis- 
tance in this region. Far from H c the magnetoresistance 
is negative again. 

The present study cannot help to answer the question 
whether there is a phase transition at H c or not be- 
cause the perturbation theory does not work very close 
to H c . Investigation of this problem can bring a new 
understanding of the granular materials. 

In conclusion, we have shown that superconducting 
fluctuations in granular metals at low temperature and 
high magnetic field lead to a suppression of the density 
of states. The well known Aslamazov-Larkin and Maki- 
Thompson corrections vanish in the limit T — > and all 
this leads to the reduction of the conductivity. 
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